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FOREWORD 
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sored  by  the  Office  of  Scientific  Research,  Mechanics  Division,  under  Contract 
AF  49(638)-952;  by  the  Office  of  Naval  Research,  Mathematical  Sciences 
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Internal  Research. 

The  material  in  this  report  was  presented  at  the  AGARD  Meeting  on 
High  Temperature  Aspects  of  Hypersonic  Fluid  Dynamics  in  Brussels,  Belgium, 
3-6  April  1962. 
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ABSTRACT 


Recent  developments  in  shock  tunnel  techniques  have  permitted  experi¬ 
mentation  in  hypersonic  airflows  having  very  high  enthalpies.  In  the  present 
paper,  the  Cornell  Aeronautical  Laboratory  six-foot  shock  tunnel  will  be  des¬ 
cribed,  the  effects  of  nonequilibrium  flow  in  such  a  device  will  be  assessed, 
and  high-enthalpy  experiments  involving  nozzle  ionization  will  be  discussed. 

The  six-foot  shock  tunnel  at  the  Cornell  Aeronautical  Laboratory 
utilizes  heated  high-pressure  hydrogen  for  a  driver  gas  to  obtain  hypersonic 
airflows  at  high  stagnation  temperatures  and  pressures.  Reflected  shock 
operation  is  used  with  the  tailored-interface  technique  to  maximize  the  dur¬ 
ation  of  steady  flow.  A  significant  portion  of  re-entry  flight  conditions  can 
be  duplicated  for  various  hypersonic  vehicles;  however,  partial  simulation 
may  be  used  to  cover  a  far  wider  range  of  velocity  and  altitude.  For  measuring 
ionization  levels  and  collision  frequencies,  microwave  diagnostics  have  been 
recently  employed. 

Owing  to  the  high  degrees  of  initial  dissociation  and  the  large  expansion 
of  the  test  gas,  it  has  been  necessary  to  study  the  influences  of  nonequilibrium 
chemistry  and  ionization  on  the  conditions  in  the  test  section.  To  this  end,  a 
general  numerical  code  has  been  developed  for  an  IBM  704  computer  in  order 
to  solve  quasi-one-dimensional  flows  with  coupled  rate  processes.  Such 
theoretical  studies  are  essential  for  planning  and  interpreting  experiments  in 
high-temperature  tunnel  flows.  Solutions  have  been  obtained  for  airflows 
with  atoms,  molecules,  electrons,  and  ions  reacting  simultaneously.  Inves¬ 
tigations  of  various  ion  «.inetic  machanisms  have  shown  that  processes  such 
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as  charge  transfer,  impact  ionization,  and  mutual  neutralization  may  play 
significant  roles  under  certain  conditions  of  hypersonic  flow  in  shock  tunnels 
and  about  blunt  bodies  in  flight. 

Preliminary  experiments  carried  out  in  the  shock  tunnel  have  dealt 
with  nonequilibrium  ionization.  Phase  shift  of  microwaves  passed  through  the 
airflow  in  the  tunnel  nozzle  has  yielded  information  on  deionization  kinetics 
and  on  the  duration  of  useful  test  flow.  Critical  frequencies  far  below  the 
microwave  band  are  indicated  for  the  gas  over  the  entire  range  of  conditions 
in  the  test  section.  These  experiments  have  indicated  that  flows  of  high 
enthalpy  can  be  generated  by  shock  tunnels  for  the  detailed  study  of  the  ion¬ 
ized  shock  layer  on  a  hypersonic  vehicle. 
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I.  INTRODUCTION 


Recent  interest  in  phenomena  associated  with  objects  entering  the  earth's 
atmosphere  at  high  velocities  has  stimulated  research  in  the  field  of  nonequili¬ 
brium  flows.  In  many  cases  a  more  detailed  knowledge  of  physico-chemical 
processes  is  needed  than  for  purely  aerodynamic  effects  because  many  of  the 
composition  changes  which  have  negligible  influence  on  the  energetics  of  the 
flow  may  be  significant  in  controlling  the  electromagnetic  properties  of  the  gas. 
This  is  true  for  the  active  and  passive  radiative  behavior  ranging  in  wavelength 
from  several  centimeters  through  the  ultraviolet  region. 

In  the  present  work,  attention  is  directed  to  the  aspects  of  the  flow  which 
determine  its  response  to  irradiation  in  the  regimes  appropriate  to  communi¬ 
cations  and  radar.  Ionization  phenomena  are  of  particular  interest  here;  spe¬ 
cifically,  the  spatial  distributions  of  electron  density  and  collision  frequency 
must  be  known  for  predicting  effects  in  the  radio-frequency  bands.  In  cases  of 
weak  signal  interactions,  the  determination  of  these  quantities  is  a  problem 
distinct  from  that  of  solving  the  equations  governing  transmission,  reflection, 
and  attenuation  of  electromagnetic  waves.  ^  Experimental  information  regard¬ 
ing  the  aerothermochemical  phenomena  is  less  extensive  than  that  pertaining 
purely  to  the  propagation  of  electromagnetic  waves  in  plasmas. 

Measurements  of  the  effects  brought  about  by  ionized  flows  on  radio  prop¬ 
agation  have  been  performed  by  means  of  flight  tests  and  by  laboratory  simu¬ 
lation  of  certain  features  characterizing  flight  conditions.  The  expense  and 
difficulty  of  flight  tests  limit  the  number  of  experiments  which  may  be  per¬ 
formed,  thereby  compromising  the  feasibility  of  exploring  fully  the  influences 
of  flight  speed,  altitude,  and  antenna  location.  As  has  been  customary  in  other 
aerodynamic  studies,  these  disadvantages  have  led  to  the  use  of  laboratory 
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test  facilities.  Clearly,  both  types  of  experimentation  are  necessary  in  order 
to  obtain  a  complete  phenomenological  picture.  Laboratory  propagation  measure¬ 
ments  have  been  reported  for  microwave  transmission  from  a  model  in  an  arc- 

2 

driven  hypervelocity  wind  tunnel.  Interpretation  of  such  data  has  been  difficult 
because  of  electrode  contamination  of  the  test  gas  and  arc  interference  effects 
on  the  instrumentation.  Another  experimental  scheme  is  the  placement  of  an 

3 

antenna  in  a  low- speed  free  jet  composed  of  seeded  combustion  gases.  Al¬ 
though  the  flow  field  and  gas  kinetic  conditions  encountered  in  flight  are  not  simu¬ 
lated  in  this  experiment,  the  ionization  levels  can  be  adjusted  by  control  of  the 

potassium  seeding.  For  electromagnetic  properties  of  wakes  behind  high- 

4  5 

velocity  projectiles,  free  flight  ranges  *  have  been  employed.  Scaling  tech¬ 
niques  such  as  those  given  by  Gibson  and  Marrone^  may  be  extended  to  aid  the 
interpretation  of  range  results. 

The  use  of  the  shock  tunnel  for  research  in  non  equilibrium  ionization  phe¬ 
nomena  will  be  discussed  in  this  paper.  For  several  years,  shock  tunnel  fa¬ 
cilities  have  been  in  use  at  the  Cornell  Aeronautical  Laboratory.  The  develop¬ 
ment  of  this  tool  for  hypersonic  research  has  been  described  in  some  detail  in 

7 

a  recent  report  by  Hertzberg  et  al.  In  the  present  work,  therefore,  only  a 
brief  description  of  the  facility  and  instrumentation  will  be  given  (Section  II). 
Following  this  descriptive  section,  the  flow  physics  of  the  nonequilibrium  phe¬ 
nomena  governing  charged  particle  distributions  will  be  discussed  (Section  III), 
The  effects  of  various  ion  kinetic  rate  processes  have  been  studied  both  for 
flows  in  a  shock  tunnel  nozzle  and  about  a  hypersonic  blunt  body.  Regimes  are 
delineated  in  which  charge  transfer  reactions  and  impact  ionization  of  molecules 
may  be  important.  Other  conditions  are  found  where  the  simple  kinetic  model 
involving  only  the  associative  ionization  (N  +  Os=2SiNO^  +  e")  is  adequate. 
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Finally,  the  results  of  recent  shock  tunnel  experiments  will  be  presented 
(Section  IV).  In  these  experiments,  microwave  diagnostic  techniques  have  been 
employed  for  measuring  electron  densities  and  electron  collision  frequencies. 
The  data  have  yielded  information  on  the  duration  of  useful  test  flow. 
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II.  SHOCK  TUNNEL  FACILITY  AND  INSTRUMENTATION 
A.  Description  of  the  Facility 

Basically,  the  shock  tunnel  is  a  short  duration  blowdown  wind  tunnel  sup¬ 
plied  by  gas  which  has  been  processed  by  shock  wave  compression.  A  nozzle 
system  fitted  to  the  downstream  end  of  a  shock  tube  serves  to  expand  the  gas  so 
that  stored  energy  is  converted  to  directed  kinetic  energy.  Either  of  two  modes 
of  operation  may  be  employed:  a)  the  direct  mode  in  which  the  incident  shock 
wave  passes  directly  through  the  nozzle  or  b)  the  reflected  mode  in  which  the 
stagnant  gas  processed  by  the  reflected  shock  wave  feeds  the  nozzle.  At  Cornell 
Aeronautical  Laboratory,  the  latter  mode  is  used  in  conjunction  with  the 

g 

"tailored-interface"  principle  in  order  to  maximize  testing  time.  Equilibrium 
interface  operation  is  used  for  obtaining  higher  enthalpy  conditions  than  would 
be  obtained  with  tailoring. 

For  the  purposes  of  duplicating  flight  conditions  in  which  ionization  phe¬ 
nomena  are  of  interest,  the  shock  tunnel  must  be  designed  to  generate  high- 
enthalpy  flows.  Additional  requirements  include  the  ability  to  control  the  im¬ 
purity  levels  in  the  test  gas  within  acceptable  limits  and  the  capability  of  ac¬ 
commodating  relatively  large-scale  models  in  the  test  section.  Impurities  must 
be  held  to  concentrations  low  enough  that  their  effect  in  producing  or  attaching 
electrons  is  negligible  compared  with  the  ionization  effects  from  the  components 
in  air.  Since  the  degrees  of  ionization  which  are  significant  in  experiments  with 
radio-frequency  wave  propagation  may  be  small,  a  high  purity  is  sometimes  a 
rather  stringent  requirement.  Large  models  are  necessary  because  complete 
similitude  can  seldom  be  attained  for  the  combination  of  aerodynamic,  nonequi¬ 
librium,  and  electromagnetic  phenomena.  In  this  connection,  the  introduction 
of  the  electromagnetic  effects  is  chiefly  responsible  for  the  difficulties  since  the 
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simultaneous  scalability  of  aerodynamic  and  nonequilibrium  effects  has  been 
theoretically  demonstrated^’*^’  for  a  wide  range  of  condition.  The  additional 
requirements  for  electromagnetic  similitude  in  hypersonic  flows  have  been  dis¬ 
cussed  in  a  previous  paper.  ^ 

In  an  effort  to  meet  the  requirements  touched  upon  in  the  preceding  para¬ 
graph,  a  six-foot  hypersonic  shock  tunnel  has  been  developed  at  Cornell  Aero¬ 
nautical  Laboratory.  Figure  1  illustrates  the  physical  layout  of  the  tunnel.  The 
driver  tube  is  equipped  with  a  heating  jacket  and  is  designed  for  the  use  of  hydro¬ 
gen  at  temperatures  up  to  750®F  at  a  pressure  of  2000  atm.  This  tube  is  14 
feet  long  and  has  an  inside  diameter  of  4. 75  inches.  The  driven  tube,  which  is 
28  feet  in  length  and  3.  00-inch  inside  diameter,  is  connected  to  a  multi-stage 
nozzle  system.  Figure  2  shows  a  detailed  view  of  the  nozzle  system  in  partial 
section.  The  two-dimensional  expansion  contour  of  the  primary  nozzle  is  de¬ 
signed  for  wave  cancellation  in  equilibrium  airflow.  The  equilibrium  Mach 
number  at  the  exit  of  this  nozzle  is  approximately  4.  5.  The  throat  area  of  the 
primary  nozzle  is  chosen  so  that  the  exhaustion  of  the  reservoir  test  gas  as 
determined  by  the  quasi- steady  mass  flow  rate  occurs  at  the  time  of  arrival  of 
the  first  expansion  wave  from  the  driver  at  the  nozzle  entrance.  A  10'  flow 
turning  section  removes  solid  particles  from  the  air  stream  by  centrifugal 
action  and  bleeds  the  boundary  layers.  The  final  stage  is  a  conical  nozzle  (10® 

-  half  angle)  in  which  the  flow  expands  to  the  six-foot  test  section  and  receiver 
vessel  shown  in  the  foreground  of  Fig.  1.  Interchangeable  liners  for  the  pri¬ 
mary  nozzle  and  entrance  lips  for  the  terminal  nozzle  provide  means  for  vary¬ 
ing  the  test  Mach  number  range.  The  conical  portion  between  the  one-foot  dia¬ 
meter  station  and  the  six-foot  test  section  is  made  of  bonded  glass  fiber  having 
a  thickness  for  optimal  transmission  of  X-band  microwave  radiation  at  9396 
Mc/sec . 
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B.  Tunnel  Performance 


Figure  3  shows  the  range  of  flight  duplication  conditions  theoretically  ac¬ 
cessible  with  the  tunnel.  The  altitude -speed  characteristics  are  particularly 
well  suited  for  the  study  of  flow  conditions  about  re-entering  satellites  and 
boost- glide  vehicles. 

The  measured  performance  of  the  shock  tube  using  heated  hydrogen  as  the 
driver  gas  and  air  as  the  test  gas  is  shown  in  Fig.  4  and  represents  the  current 
high-performance  limit.  This  figure  shows  a  plot  of  stagnation  enthalpy  and 
reflected- to-driver  pressure  ratios  versus  shock  Mach  number.  The  thermo¬ 
dynamic  state  of  the  gas  in  this  reflected  shock  region  is  determined  from 
measurements  of  the  incident  shock  velocity  and  the  stagnation  pressure.  It  has 
been  found  that  for  shock  tube  operation  near  the  tailored-operating  point  that 
the  waves  generated  due  to  the  mismatch  are  quickly  attenuated  so  that  a  steady 
reservoir  condition  is  in  fact  attained.  As  the  waves  generated  are  very  weak, 
the  final  state  is  determined  by  assuming  an  isentropic  process  to  occur  from 
the  region  behind  the  first  reflected  shock  to  the  final  steady  state  condition. 

In  order  to  obtain  a  sufficient  duration  of  uniform  reservoir  conditions,  the 
shock  tube  is  operated  near  the  tailored-interface  conditions.  Contamination  of 
the  test  gas  by  mixing  with  the  driver  gas  is  minimized  by  operation  under 

s 

conditions  where  the  interface  separating  the  driver  and  driven  gases  is  stable. 
In  this  context,  stability  refers  to  that  condition  where  the  density  of  the  hydro¬ 
gen  gas  is  less  than  that  of  the  shocked  air  at  the  interface  zone  so  that  any 
mixing  due  to  instabilities  does  not  contaminate  the  region  of  the  nozzle  inlet. 

Calibration  studies  of  the  nozzle-test  section  system  have  included  pitot 
pressure  surveys  across  the  tunnel  using  a  multiple-probe  rake.  The  pitot 
pressure  distributions  at  the  test  section  for  a  range  of  reservoir  conditions 
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are  shown  in  Fig.  5.  These  data  were  taken  with  the  2-5/8"  inlet  diameter  of 
the  second- stage  nozzle. 

Information  on  flow  angularity  with  respect  to  the  10°  nozzle  axis  has  been 
obtained  from  the  results  of  heat  transfer  studies  on  a  slender  cone^^  following 
techniques  described  in  Ref.  12.  The  results  of  these  tests  indicate  the  center- 
line  flow  at  the  test  section  station  to  have  a  total  angularity  of  less  than  1*. 

C.  In  s  tr  um  entat  i  on 

Methods  for  the  accurate  collection  of  data  under  conditions  of  high  ve¬ 
locity,  high  temperature,  and  short  test  flow  duration  are  necessary  in  shock 
tunnel  research.  Cathode-ray  oscilloscopes  provide  adequate  means  of  record¬ 
ing  pressure  and  heat  transfer  information.  Rapid  response  instrumentation 
for  measuring  the  wide  range  of  pressures  encountered  in  shock  tunnel  research 
includes  modified  microphones  and  a  variety  of  piezoelectric  crystal  trans- 
ducers.  The  design  of  the  transducer  porting  incorporates  shielding  from  the 
thermal  radiation  of  the  gas  flow  in  such  a  way  that  excessive  blockage  does  not 
seriously  increase  the  over-all  response  time.  Transient  surface  temperatures 
are  measured  with  thin-film  resistance  thermometers  mounted  on  dielectric 
backing  material.  Specifically,  the  gauges  consist  of  platinum  films  of  the 
order  of  1000  X  thick  painted  or  sputtered  on  a  Pyrex  backing  block.  Heat  fluxes 

2  4 

are  recorded  directly  from  an  analogue  network  which  receives  the  tempera- 

15  . 

ture  signal  from  the  thermometer.  A  pulse  method  is  employed  for  determin¬ 
ing  the  thermal  properties  of  the  backing  material.  Thin  insulative  coatings  are 
applied  to  the  outer  surface  of  the  thin  film  in  order  to  block  any  electrical  con¬ 
duction  paths  through  the  ionized  gas  flowing  over  the  surface.  Detailed  dis¬ 
cussions  of  the  theory  and  application  of  the  aerodynamic  instrtimentation  may 
be  found  in  Ref.  7. 
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The  thin-film  gauges  have  been  found  to  respond  to  microwave  radiation 
as  well  as  aerodynamic  heat  transfer.  In  order  to  separate  the  two  effects, 
gauges  are  mounted  both  on  the  exterior  and  interior  surfaces  of  the  model,  and 
the  wall  of  the  model  is  designed  for  microwave  transparency.  Glick^^  has  in¬ 
vestigated  the  microwave  response  of  these  gauges  theoretically  and  experi¬ 
mentally.  A  heat  transfer  analysis  was  carried  out  for  a  two-layer  slab;  the 
first  layer  is  assumed  to  be  bounded  on  one  side  by  an  adiabatic  surface  and  on 
the  other  by  the  interface  with  the  second  layer  which  is  assumed  to  be  semi¬ 
infinite.  Uniform  volumetric  heat  generation  is  specified  for  the  first  layer. 
Solutions  of  the  heat  equations  and  substitution  of  the  properties  of  the  thin 
platinum  film  and  the  Pyrex  backing  show  that  the  temperature  in  the  film  is  es¬ 
sentially  uniform  for  times  greater  than  lO’^  seconds.  For  these  times,  the 
temperature  time  variation  of  the  film  approaches  that  at  the  surface  of  a  semi- 
infinite  homogeneous  solid  with  constant  heat  input  equal  to  the  energy  generation 
in  the  gauge  per  unit  area  of  film  surface.  The  antenna  characteristics  of  the 
gauge  and  its  equivalent  circuitry  were  also  analyzed  in  Ref.  17.  It  was  found 
that  the  generation  term  was  directly  proportional  to  the  electromagnetic  energy 
flux.  The  proportionality  factor  can  be  determined  theoretically  on  the  basis  of 
the  analytical  results.  An  experiment  was  performed  to  determine  the  depend¬ 
ence  of  equivalent  heat  flux  on  the  angle  between  the  electric  field  direction  and 
the  dipole  direction  characterizing  the  gauge  if  it  were  considered  to  be  an  anten¬ 
na.  The  results  verified  the  fact  that  the  gauge  acts  as  a  dipole  antenna  rather 
than  a  simple  absorber,  which  would  not  exhibit  any  angular  dependence.  In 

.  7 

Ref.  17  an  example  is  given  for  a  platinum  thin-film  bolometer  of  thickness  10 

_3 

meters  and  length  of  about  3x10  meters.  The  useful  range  of  this  gauge  as 

2 

a  bolometer  is  estimated  to  be  0,0Z5  to  23  watts/cm  of  incident  microwave 
radiation  flux. 
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The  electronic  components  used  for  both  aerodynamic  and  microwave  in¬ 
strumentation  are  standard  for  the  most  part,  and  will  be  described  diagram- 
matically  in  the  context  of  the  experimental  discussion. 


m.  THEORETICAL  STUDIES  OF  NONEQUILIBRIUM  PHENOMENA 


A.  Delineation  of  Problems 

Some  of  the  most  important  features  in  the  interpretation  of  data  obtained 
in  high- enthalpy  shock  tunnel  research  involve  nonequilibrium  phenomena.  Be¬ 
cause  of  the  large  degrees  of  dissociation  and  ionization  initially  existing  in  the 
test  gas  before  expansion,  the  aerothermochemistry  of  the  nozzle  flow  must  be 
analyzed  in  order  to  aetermine  conditions  in  the  test  section  of  the  tunnel.  Fur¬ 
thermore,  certain  aspects  of  the  hypersonic  flow  about  a  model  (or  a  flight  ve¬ 
hicle)  are  dominated  by  nonequilibrium  effects.  In  the  present  work  the  influ¬ 
ence  of  coupled  rate  processes  on  the  distributions  of  charged  particles  will  be 
emphasized  because  the  chemical  aspects  of  the  problem  have  been  discussed 
in  detail  in  previous  papers  (e.  g. ,  Ref.  10).  Kinetic  mechanisms  of  varying 
degrees  of  complexity  will  be  examined  in  Subsection  B  in  order  to  evaluate  ex¬ 
isting  simplifications  and  their  effects  on  electron  distributions  in  particular. 
After  the  computational  models  are  established,  results  from  the  various  models 
will  be  compared  in  Subsection  C  for  nozzle  flows  at  conditions  appropriate  to 
shock  tunnel  testing.  In  Subsection  D,  the  ionization  histories  will  be  examined 
for  the  stagnation  streamtube  in  a  hypersonic  blunt-body  flow  for  various  kinetic 
models. 

B .  Basis  of  Ionization  Kinetic  Models 

In  most  of  the  hypersonic  flow  conditions  considered  here,  the  rate  pro¬ 
cesses  involving  the  neutral  particles  dominates  the  thermodynamic  history  of 
a  fluid  element  and  the  ion  kinetics  are  governed  by  this  history.  The  neutral 
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particle  reaction  model  employed  in  Ref.  10  is  used  for  all  cases  (see  Table  I). 
The  use  of  this  set  of  reactions  is  fully  discussed  in  Ref.  10. 

A  major  portion  of  previous  work  on  nonequilibrium  ionization  in  air  has 
been  limited  to  associative  ionization  (and  the  reverse,  dissociative  recombi¬ 
nation)  types  of  reactions  because  these  are  energetically  favored.  Exceptions 
to  this  include  Refs.  18  and  19.  In  Ref.  18  two  charge  transfer  reactions  were 
added  in  shock  wave  calculations  and  in  Ref.  19  the  added  effects  of  atomic  ions 
governed  only  by  three -body  Thomson  recombination  were  examined  in  expansion 
calculations.  Both  additional  effects  are  studied  simultaneously  in  this  paper. 
For  shock  tunnel  expansions  from  relatively  high-pressure  reservoir  conditions, 
the  negative  oxygen  atomic  ion  is  added  in  one  case  to  determine  its  effect  on 
teat  flow  electron  concentrations.  When  O  is  included,  collisional  detachment 
and  mutual  neutralization  reactions  are  included.  Radiative  and  dielectronic 
recombination  processes  are  omitted  for  the  reasons  given  in  Ref.  19;  i.e.  , 
faster  competitive  recombination  paths  exist  for  the  range  of  conditions  chosen. 
Table  2  summarizes  the  ionization  reactions  and  the  corresponding  rate 

constants  which  have  been  employed  in  the  present  studies.  Reactions  10-12 

-I-  +  4- 

are  dissociative  recombination  of  NO  ,  N^  and  O  respectively  with  elec- 
trons.  Reactions  13-22,  40-49  are  three-body  recombinations,  reactions  23- 
33  are  charge  transfers,  reactions  34-38  are  mutual  neutralizations,  and  re¬ 
action  39  is  three-body  attachment.  The  rate  constant  used  for  reaction  10 
is  that  obtained  experimentally  by  Lin^^  and  reported  in  Ref.  18.  Axial  micro- 
wave  reflection  experiments^^  carried  out  in  a  shock  tube  at  CAL  yield  values 
for  ionization  relaxation  distance  of  the  same  magnitude  as  that  predicted  from 

rate  data  given  in  the  earlier  work  cited.  Theoretical  calculations  of  the  cross 

22 

section  have  also  been  reported  recently.  The  value  obtained  for  1  ev. 
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electrons  is  an  order  of  magnitude  greater  than  that  given  in  Ref.  18.  It  should 
be  pointed  out  that  the  theoretical  estimate  is  for  monoenergetic  electrons,  and 
is,  therefore,  not  specifically  comparable  with  experimental  results  for  ther- 
malized  electrons.  Microwave  phase- shift  measurements  of  ionization  in  a 
nozzle  flow  are  being  employed  at  the  time  of  this  writing  to  determine  the  rate 
constant  for  reaction  10.  Early  indications  are  that  is  the  same  order  of  mag¬ 
nitude  as  the  value  listed  in  Table  2.  This  shock  tunnel  nozzle  experiment  and 
its  theoretical  interpretation  are  the  subjects  of  Section  IV  of  this  paper.  Rate 
constants  for  reactions  11  and  12  are  baaed  on  theoretical  estimates  given  in 
Ref.  18.  Reactions  13-22  and  40-49  are  three-body  Thomson  recombination 
reactions.  The  rate  data  for  these  reactions  were  computed  theoretically  ac¬ 
cording  to  the  classical  Thomson  theory  modified  where  necessary  for  inelastic 
collisions  of  electrons  with  molecules.  The  inelasticity  factors  were  taken  to 
be  constant  and  were  obtained  by  averaging  over  the  energy  range  of  interest. 

Electron-ion  recombination  reactions  (10-22)  constitute  parallel  paths  for 
electron  removal  which  are  coupled  only  through  the  electron  concentration. 

Over  most  of  the  density  range  of  interest  in  hypersonic  flight  and  in  shock 
tunnel  operation,  the  binary  reactions  (10-12)  are  by  far  the  most  active  paths 
in  electron  removal  and  their  reverse  reactions  likewise  are  the  most  signifi¬ 
cant  in  electron  production.  The  atomic  ions  may  be  neutralized  relatively 
efficiently  by  means  of  charge  transfer  processes  of  the  types 

A"*"  +  BC  AB'*’  +  C 
+  BC  +  BC'*’ 

which  are  required  for  completion  of  the  coupling  through  atomic  ions.  A,  B, 
and  C  are  symbols,  any  of  which  may  be  N  or  O,  These  classes  of  reaction 
produce  diatomic  ions,  which  have  a  much  higher  probability  of  reaction  with 
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electrons  than  the  atomic  ions  because  of  the  dissociative  processes.  As  brought 
out  in  the  classical  chain  reaction  theory  of  chemical  kinetics,  certain  chain 
carriers  may  have  a  strong  influence  on  over-all  kinetics  despite  their  relatively 
small  equilibrium,  or  quasi- stationary  concentration.  For  this  reason,  the  ef¬ 
fect  of  charge  exchange  reactions  on  electron  densities  v/ill  be  considered  in  the 
« 

present  work.  Although  experimental  rate  data  are  unavailable  for  the  charge 
transfer  reactions  in  the  temperature  range  of  interest  here,  conservative  esti¬ 
mates  have  been  employed.  The  rate  constants  for  reactions  23-33  are  those 

given  in  Ref.  18  which  are  based  on  the  assumption  of  cross  sections  of  order 

a  1/2 

A'3.0  for  all  charge  transfer  reactions  (hence  the  T  '  dependence).  8.^  is 

1 8 

the  first  Bohr  radius.  Such  an  assumption  is  reported  to  agree  with  unpub¬ 
lished  experimental  data  for  reaction  25. 

23 

The  theory  of  Eyring,  Hirschfelder ,  and  Taylor  gives  a  reaction  cross 
section  inversely  proportional  to  the  relative  speed  of  collision  partners  so  that 
the  rate  constant  would  be  independent  of  temperature.  This  theory  was  devel¬ 
oped  for  particles  which  obey  a  Maxwellian  velocity  distribution  and  which  have 
no  permanent  dipole  moment.  The  energy  of  activation  is  considered  to  be  that 
at  the  maximum  in  a  surface  defined  by  the  difference  between  the  rotational 
energy  of  the  colliders  and  the  polarization  energy  of  the  encounter.  The  rate 
constant  is  determined  to  a  constant  (a  transmission  coefficient)  by  carrying 
out  the  collision  integral  to  determine  the  rate  of  passage  over  the  barrier  in 
the  energy  surface.  For  a  transmission  coefficient  of  unity,  the  expression  for 
the  rate  constant  is  as  follows 


where  oc  is  the  polarizability,  e  the  electonic  charge  and^  the  reduced 
mass  of  the  colliding  pair.  Table  3  summarizes  the  polarizability  data  used 
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and  the  rate  constants  obtained  with  the  theory  of  Ref.  23.  Reactions  26,  27,  31, 
and  32  are  omitted  because  nitric  oxide  has  an  electric  dipole  moment,  ,  of 
1.6  X  10  esu,  and  the  theory  of  Ref.  23  is  not  strictly  applicable.  An  estimate 
for  the  polarizability  may  be  obtained  by  using  the  formula  developed  in  Ref.  24, 
oc  =  oCj  +  with  as  determined  in  Ref.  25  (1.68  x  10  cm^). 

At  6000 *K  the  permanent  moment  contribution  given  by  this  formula  is  of  the 
order  of  1%  of  the  induced  moment  because  of  thermal  randomization.  Hence, 
oC  zi  oCg  ,  and  the  charge  transfer  rate  constants  for  reactions  26,  27,  31,  and 
32  may  be  expected  to  be  of  the  same  order  as  the  others  in  the  limit  of  this 
approximation. 

Numerical  estimates  of  the  charge  transfer  rate  constants  have  also  been 

26  27  28 

given  by  Beran  for  5000  to  6000 “K  and  by  Meyerott,  Krasovskii,  and 

29 

Nawrocki  for  upper  atmospheric  conditions  (presumably  several  hundred  *K). 

The  various  estimates  based  on  theory  and  observation  are  spotted  on  a  graph  of 

rate  constant  vs  temperature  in  Fig.  6.  Note  that  the  line  based  on  the 

cross  section  is  conservative  compared  with  all  other  estimates  except  that  of 

+  29 

Ref.  29.  Nawrocki  gives  rather  wide  uncertainty  margins  (log^^k  =111  1). 
Therefore,  this  source  may  be  weighted  less  than  the  others.  The  use  of  con¬ 
servative  estimates  of  k22  33  will  minimize  the  effects  which  they  have  on  the 
system.  This  procedure  is  necessary  because  of  the  uncertainties  in  the  rate 
constants. 

In  the  density  regimes  encountered  at  the  early  portion  of  nozzle  expan¬ 
sions  typical  of  shock  tunnel  testing,  significant  amounts  of  O  might  be  ex¬ 
pected  because  oxygen  atoms  are  abundant  and  have  relatively  high  electron 
affinity  (1.46  ev).  Either  the  affinities  of  the  other  species  present  are  much 
smaller  than  that  of  the  oxygen  atom  or  the  neutral  species  in  question  are 
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present  in  much  lower  concentration.  Therefore,  all  other  negative  ions  are 
neglected  in  this  work.  The  reactions  involving  this  species  are  Nos.  34-39  ih 
Table  Z.  All  but  the  last  reaction  of  this  group  are  mutual  neutralizations. 

Since  the  interaction  is  governed  by  the  long-range  coulombic  potential  and  the 
excess  energy  may  be  carried  off  by  recoil  (as  well  as  excitation),  mutual  neu¬ 
tralization  is  an  extremely  efficient  type  of  reaction.  Reaction  39  is  a  three- 
body  attachment  where  M  represents  any  third  body.  This  reaction  is  a  close 
analogue  to  three -body  neutral  atom  recombination  -  in  fact  its  rate  constant  at 
300*K  is  of  the  same  order  of  magnitude  as  that  for  the  corresponding  reaction 
involving  neutral  species.  Dissociative  attachment  to  and  NO  are  omitted 
because  these  reactions  have  moderate  activation  energies.  Hence,  their  rates 
would  vanish  exponentially  with  decreasing  temperature.  Rate  constants  sug¬ 
gested  in  Ref.  29  were  employed  for  reactions  34,  36,  and  39.  The  negative 
one-half  power  temperature  dependence  for  reactions  34-38  (inclusive)  was  also 

given  in  Ref.  29.  Rate  data  for  reaction  35  were  obtained  from  Massey's 
30 

article  on  gaseous  ion  recombination,  and  those  for  reactions  37  and  38,  from 
Mitra..^^ 

For  ease  in  future  reference,  the  various  reaction  kinetic  models  will  be 
numbered  as  follows: 
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Model 

Number 

Species 

Reactions* 

Remarks 

1 

1  -  8 

1  -  10 

Basic  model 
(NO"'’  +  e"  +  0) 

2 

1  -  12 

1  -  22 

Add  three-body  recombi¬ 
nation  of  atomic  ions 

3 

1  -  12 

1  -  32 

Add  charge  transfer 

4 

1  -  12 

1  -  33 

Add  charge  transfer 

5 

1  -  13 

1  -  39 

Add  mutual  neutralization 
and  attachment 

6 

1  -  12 

1  -  22,  40  -  49 

40  -  49  are  three-body  re¬ 
combination  involving  mo¬ 
lecular  ions 

7 

1  -  12 

1  -  12,  23  -  33,  40  -  49 
15  -  17,  20  -  22 

C,  Shock  Tunnel  Nozale  Flows 

At  the  shock  tunnel  condition  considered  here,  (T^  s  SOOO’K,  =  100  atm), 
the  variation  of  aerotherxnodynamic  quantities  (p, ^  ,T  ,  u  *  and  M)  are  princi¬ 
pally  controlled  by  the  interaction  of  the  coupled  chemistry  of  the  neutral  species 
with  the  dynamics  of  the  expansion  process.  The  method  of  computation  and  the 
first  results  for  hypersonic  nozzle  airflow  governed  by  coupled  chemical  re¬ 
actions  appear  in  Ref.  32.  Further  analysis  and  discussion  of  these  results  in 
the  context  of  hypersonic  tunnel  testing  appear  in  Ref.  10.  In  the  rapid  expan¬ 
sion  of  the  high-enthalpy  flow  through  a  shock  tunnel  nozzle,  various  chemical 
degrees  of  freedom  may  exhibit  nonequilibrium  lags  because  molecxilar  collisions 
become  lees  and  less  frequent.  Although  vibrational  lag  processes  are  believed 

*  Species  and  reaction  numbers  are  given  in  Table  2. 
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to  have  little  effect  on  the  gasdynamical  history,  the  freezing  of  various 
chemical  processes  can  withhold  rather  large  fractions  of  energy  from  the  flow. 
Rotational  and  translational  relaxation,  being  relatively  efficient  processes,  can 
safely  be  considered  to  equilibrate  throughout  the  entire  expansion  history.  The 
fraction  of  oxygen  atoms  tends  to  become  constant  rather  early  because  the 
ternary  recombination  reactions  freeze  first  and  no  efficient  binary  paths  exist 
for  removing  oxygen  atoms.  On  the  other  hand,  nitrogen  dissociation  is  sub¬ 
stantially  equilibrated  when  sufficient  nitric  oxide  is  present  to  participate  in 
the  fast  binary  reactions  (8  and  9)  which  remain  active  after  the  ternary  re¬ 
actions  (1  -  7B)  freeze. 

Nitric  oxide  molecules  undershoot  the  equilibrium  curve  when  insufficient 
is  present  to  replenish  NO  through  reaction  8  as  NO  is  consumed  by  re¬ 
action  9.  If  nitric  oxide  is  seriously  depleted,  the  nitrogen  atom  fraction 
suddenly  halts  its  descent.  Otherwise,  it  levels  out  gradually  as  the  two-body 

reactions  become  ineffectual.  The  chemical  composition  histories  for  the  shock 

32 

tunnel  cases  considered  for  the  ionization  studies  is  shown  in  Fig.  7.  The 
nitrogen  atom  concentration  drops  until  the  exchange  reactions  freeze  in  the 
region  of  NO-undershoot. 

As  shown  in  Ref.  32,  collision  frequency  of  N-atoms  with  and  NO  is 

sufficiently  high  at  p^  =  1000  atm  to  equilibrate  nitrogen  dissociation  throughout. 

32 

Table  4  summarizes  the  gasdynamic  consequences  of  chemical  nonequilibrium. 
For  more  complete  discussion  of  the  behavior  of  neutral  particles,  the  reader 
should  consult  Refs.  10  and  32. 

Previous  work  on  nonequilibrium  ionization  in  expanding  air  flows  appears 
in  Refs.  19  and  33.  Ion  kinetic  model  1  with  the  addition  of  Thomson  recombi- 
nation  of  (e  +  NO  )  is  the  basis  of  approximate  freezing  solutions  utilizing 
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a  freezing  criterion.  Electromagnetic  wave  penetration  depths  in  the  test  flow 
of  the  CAL  six-foot  tunnel  are  also  shown  in  Ref.  33,  Two  conclusions  are 
apparent  from  these  results:  (1)  Thomson  three-body  recombination  is  com¬ 
parable  with  dissociative  recombination  in  air  only  at  the  highest  reservoir 
density  considered  (  p  ! ^5)  and,  (2)  the  flow  in  the  test  section  of  the 
tunnel  should  be  expected  to  be  transparent  to  UHF  and  microwave  radiation. 

As  will  be  described  in  a  later  section,  this  conclusion  has  been  borne  out 
experimentally. 

The  more  recent  work  on  ionization  nonequilibrium  in  airflows  is  covered 
in  Ref.  19.  Exact  solutions  are  compared  with  approximate  solutions  and  it  is 
shown  that  for  certain  flow  conditions,  the  electron  fraction  does  not  approach 
a  constant  value  asymptotically,  but  rather  decreases  as  some  power  of  dis¬ 
tance  dependent  upon  frozen  chemical  composition  and  nozzle  geometry,  For 
completely  frozen  ionization,  the  d.c.  conductivity  reversal  is  demonstrated 
for  pure  and  for  cesium- seeded  airflows.  Ion  kinetic  model  1  was  employed 
for  the  nozzle  flow  solutions  in  Ref,  19. 

In  the  present  work,  ion  kinetic  models  1,  3  and  5  will  be  compared  for 

the  reservoir  state  of  (T  =  8000®K,  p  =100  atm).  Model  1  consists  only  of 

o  o 

the  fastest  electron  removal  mechanism  because  electrons  are  completely 
counterbalanced  by  diatomic  ions  capable  of  reacting  in  two-body  recombinations, 
whereas  in  model  3,  the  electrons  have  positively  charged  counterparts  made 
up  of  both  atomic  and  molecular  ions.  If  the  removal  of  atomic  ions  were  solely 
dependent  on  reactions  with  electrons  (such  as  in  model  2),  the  ionization  system 
would  freeze  suddenly  as  soon  as  the  supply  of  molecular  ions  was  depleted  by 
dissociative  recombinations.  This  behavior  would  bear  a  close  resemblance  to 
that  of  nitrogen  dissociation  freezing,  which  occurs  in  airflows  when  the  number 
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of  nitric  oxide  molecules  diminishes  sharply.  In  model  3,  however,  reaction 
paths  (23-33)  other  than  three-body  recombination  (13-22)  are  allowed  for  the 
removal  of  the  atomic  ions. 

In  reactions  22-33,  atomic  ions  may  be  neutralized  by  acquiring  an  elec¬ 
tron  from  some  previously  neutral  collision  partner.  If  the  collision  partner  is 
diatomic,  the  positive  ion  formed  may  then  participate  in  binary  recombination 
with  electrons.  The  supply  of  diatomic  ions  depends  on  a  balance  set  up  be¬ 
tween  production  by  charge  transfer  and  consumption  by  recombination.  Thus, 
molecular  ions  are  truly  chain  carriers  in  the  neutralization  of  atomic  ions  when 
charge  transfer  occurs.  After  the  rates  of  endothermic  reactions  become  rela¬ 
tively  small,  the  chains  are  of  the  forms: 

+  BC  — >  C  +  AB'*’ 

AB'*’  +  e - -  A  +  B 

and 

A"^  +  BC  - S'  A  +  BC"^ 

BC'^  +  e' - S'  B  +  C 

These  chains  remain  effective  after  the  Thomson  processes 

M  +  A"^  +  e”  — >  A  +  M 

has  effectively  vanished  in  the  expanding  flow. 

The  detailed  chain  structure  may  be  represented  by  the  parallel  heirarchy 
of  exothermic  reactions  shown  in  Fig.  8.  From  an  energetic  point  of  view, 
seems  to  be  an  exception  to  this  rule.  Reactions  23  and  28  are  exothermic  to 
the  right  by  2. 0  and  1.0  e,  v.  respectively.  The  forward  activation  energies  are 
zero;  therefore,  the  reverse  reactions  freeze  more  rapidly  in  an  expanding  flow 
because  of  the  exponential  dependence  of  the  reverse  rate  on  temperature.  In 
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departures  from  ionization  equilibritim  then,  the  net  result  of  reactions  E3  and  28 
is  to  oppose  the  chain  for  removing  atomic  ions,  and  the  atomic  ions  themselves 
are  chain  carriers  in  a  series  of  reactions  for  ^2^  removal  in  competition  with 
reaction  12. 

With  the  additional  species  and  reactions  in  model  3  as  compared  with 
model  1,  two  alternative  questions  are  answered  by  carrying  out  exact  solutions: 
1.  Are  charge  transfer  reactions  (24-27;  29-32)  rapid  enough  to  impart  to  the 
ionization  system  the  gross  behavior  of  binary  recombination,  or  2.  Are  charge 
transfer  reactions  23  and  28  sufficiently  influential  to  oppose  the  electron  re¬ 
moval  chain  by  producing  0'*’and  at  a  rate  exceeding  removal  by  reactions 
(13-22;  24-27;  29-32). 

Figure  9  shows  the  ionization  history  in  the  nonequilibrium  nozzle  flow 
considered,  Kinetic  models  1,  3,  and  S  have  been  compared.  The  electron 
concentration  histories  for  models  1  and  3  practically  coincide  (particularly  at 
large  area  ratios).  This  shows  that  at  an  early  stage  the  charge  transfer  re¬ 
actions  effectively  eliminate  all  atomic  ions  which  initially  counterbalance  4.  3% 
of  the  electron  population.  If  these  fast  chain  reactions  were  not  considered,  the 
atomic  ions  would  recombine  only  by  the  paths  set  out  in  reactions  13-22  which 
freeze  relatively  earlier  because  of  higher  order  density  dependence.  Model  5 
contains  0~  as  an  additional  species  along  with  six  reactions.  At  the  stagnant 
reservoir  conditions,  negative  oxygen  ions  are  about  1 /20th  as  numerous  as 
electrons.  During  the  expansion,  however,  the  mutual  neutralization  reactions 
(34-38)  prevent  freezing  of  O”  which  would  occur  rather  early  if  reverse  re¬ 
action  39  (collisional  detachment)  were  the  only  removal  process  in  operation. 
This  ratio  is  maintained  or  slightly  reduced  in  the  expansion.  Compared  with 
the  case  without  O",  the  electron  population  is  reduced  by  less  than  5%. 
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From  these  comparisons  it  may  be  concluded  that  the  simple  ionization 
kinetics  represented  by  model  1  are  adequate  for  predicting  ionization  levels  in 
shock  tunnel  test  flows  for  stagnation  temperatures  up  to  8000‘’K  and  at  pres¬ 
sures  where  nonequilibrium  effects  are  significant.  Since  conservative  esti¬ 
mates  have  been  employed  for  the  charge  transfer  cross-section,  this  conclu¬ 
sion  may  well  hold  above  the  temperature  considered  here.  At  temperatures 
lower  than  8000®K,  positive  ions  other  than  NO^  are  not  significant.  Clearly, 
at  some  high  temperature,  molecular  ions  will  not  be  numerous  enough  to  main¬ 
tain  any  of  the  fast  chains  involving  charge  transfer  or  dissociative  recombi¬ 
nation,  and  freezing  of  the  ionization  may  occur  very  early. 

D.  Inviscid  Hypersonic  Flow  About  a  Blunt  Body 

The  significant  chemical  features  of  these  flows  were  discussed  in  Ref.  10 
principally  for  processes  controlling  the  distribution  of  energy  among  the  vari¬ 
ous  degrees  of  freedom.  Since  the  present  work  is  mainly  concerned  with  the 
plasma  kinetics  of  the  flow,  only  a  brief  summary  of  the  neutral  particle  aero- 
thermochemistry  will  be  given  here.  Clearly,  viscous  effects  will  play  an  im¬ 
portant  role  in  the  flow  near  the  body  surface.  However,  the  purpose  of  the 
present  work  is  the  assessment  of  ionization  kinetic  mechanisms  and  their  ef¬ 
fects  on  the  plasma  properties.  For  this  purpose,  inviscid  solutions  are  ade¬ 
quate  for  providing  the  thermodynamic  environment  typifying  high  altitude  hyper¬ 
sonic  flight. 

1.  Review  of  Neutral  Particle  Effects 

The  initial  effect  on  a  fluid  element  as  it  is  processed  by  the  detached 
shock  wave  preceding  a  blunt  body  is  rapid  excitation  of  translational  degrees  of 
freedom.  Practically  coincident  with  this  process  is  rotational  excitation.  At 
many  flight  conditions  typifying  re-entry  trajectories,  vibrational  and 
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disaociational  relaxation  are  interdependent  through  a  complex  coupling  mecha¬ 
nism.^®  As  the  molecules  in  air  dissociate  to  form  atoms,  nitric  oxide  is 
formed  by  rapid  atom  exchange  reactions.  These  reactions  approach  quasi¬ 
equilibrium  so  rapidly  that  the  nitric  oxide  concentration  overshoots  infinite 
rate  equilibrium  values;  this  occurs  to  a  great  degree  along  the  higher  entropy 
stream  surfaces  emanating  from  the  central  portion  of  the  shock  wave.  As  the 
gas  expands  and  accelerates  iu  the  flow  about  the  body,  the  temperature  and 
density  decrease,  driving  downward  the  local  equilibrivixn  atom-concentrations. 
The  presence  of  nitric  oxide  in  significant  concentrations  favors  the  removal  of 
nitrogen  atoms  through  rapid,  two-body  paths;  on  the  other  hand,  oxygen  atomic 
recombination  freeses  early.  When  the  reaction  of  N  with  NO  disilpates  the 
supply  of  nitric  oxide  more  rapidly  than  the  N  O2  reaction  (8)  replenishes  it, 
nitrogen  atom  recombination  freezes  and  nitric  oxide  undershoots  equilibrium. 
Except  for  altitudes  of  100  Kft  and  below,  nitric  oxide  becomes  deplsted  so 
rapidly  that  the  nitrogen-atom  fraction  changes  very  little  from  its  stagnation 
level  in  the  flow  near  the  surface  of  typical  re-entry  vehicles.  In  outer  Stream¬ 
lines,  freezing  may  occur  before  dissociation  is  complete  because  of  the  high 
velocity  and  rapid  expansion  Immediately  behind  the  shock. 

These  features  of  nonequilibrium  airflows  about  blunt  hypersonic  bodies 
have  been  brought  out  in  numerical  solutions  obtained  by  means  of  computing 
machine  codes.  The  roles  of  the  endothermic  processes  in  the  shock  wave  de¬ 
tachment  region  about  the  nose  cap  have  been  investigated  by  the  inverse  method 

34 

which  has  been  extended  to  coupled  nonequilibrium  flows.  The  results  of  this 

35 

work  were  described  in  a  brief  paper.  In  addition  to  the  series  of  reports 
under  Ref.  34,  detailed  results  and  discussions  of  nonequilibrium  behavior  are 
found  in  Ref.  10.  The  expansion  freezing  behavior  in  the  streamtubes  originating 
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from  the  stagnation  equilibrium  region  has  been  studied  by  means  of  an 
adaptation  ’  of  the  nozzle  flow  code  to  prescribed  pressure  distributions. 

In  the  present  work,  the  ionization  kinetics  are  studied  for  one  flight 
condition  typical  of  a  ballistic  re-entry  vehicle  (speed  =  23  Kft/sec,  altitude  = 
150  Kft)  and  one  flight  condition  typical  of  a  lunar  probe  re-entry  vehicle  (speed 
36  Kft/ sec,  altitude  =  200  Kft).  For  illustrative  purposes,  flow  field  solutions 
are  limited  to  a  stagnation  streamline  starting  from  the  detached  shock,  going 
centrally  inward  to  the  stagnation  point,  and  proceeding  around  the  body  at  the 
surface.  The  body  chosen  is  a  hemisphere- cylinder  with  a  1-foot  nose  radius 
flying  at  zero  yaw.  In  the  ballistic  re-entry  condition  the  energy  residing  in 
ionization  is  small  compared  with  the  total  enthalpy;  however,  in  the  lunar  probe 
re-entry  condition  nearly  18%  of  the  energy  is  in  ionization.  For  the  decelera¬ 
ting  portion  of  the  streamtube  flow  between  the  shock  wave  and  the  stagnation 
point,  the  normal  shock  wave  numerical  code^^  was  used  with  a  streamwise 
coordinate  transformation^  to  adapt  the  parallel  flow  case  to  the  stagnation 

point  flow.  For  the  accelerating  flow  from  the  stagnation  point  around  the  nose 

32 

and  along  the  afterbody,  an  adaptation  of  the  nozzle  flow  numerical  code  was 
used.  Modified  Newtonian  and  blast  wave  pressure  distributions  were  assumed 
for  the  nose  cap  and  afterbody  respectively  as  described  in  Ref.  10. 

2.  Ballistic  Re-eniry  Case 

Three  cases  ar-  compared  for  the  ion:zation  in  the  flow  about  a  ballistic 
re-entry  vehicle. 

Case  A  consists  of  reaction  kinetic  model  number  1  for  the  entire  flow. 

This  is  included  for  comparative  purposes  because  much  of  the  previous 

literature  on  air  ionization  is  limited  to  sijich  a  mechanism. 


Cases  A  and  B  were  previously  .shown  in  Ref.  19. 


Case  B  is  based  on  model  6  for  the  decelerating  flow  and  on  model  2  for 
the  accelerating  flow.  In  these  models,  all  of  the  singly  positive  ions  are 
included  as  well  as  electrons.  Model  6  differs  from  model  2  in  that  it 
allows  impact  ionization  of  diatomic  species.  The  reverse  reaction,  three- 
body  recombination  of  electrons  with  positive  diatomic  ions,  is  overshadowed 
by  dissociative  recombination  in  the  exothermically  lagging  condition  in  ex¬ 
panding  flow;  hence,  the  choice  of  model  2  for  that  flow.  Early  in  the  history 
of  the  shock  wave  relaxation  zone,  atoms  are  not  yet  produced  in  large 
quantities  and  translational  energy  is  large.  Although  the  probability  of 
dissociation  is  much  higher  than  that  for  ionization,  it  is  interesting  to  com¬ 
pare  the  impact  ionization  rate  with  the  associative  ionization  rate  in  the 
early  stages. 

Case  C  includes  the  addition  of  the  charge  transfer  reactions  to  Case  B. 
Model  7  is  employed  for  the  decelerating  flow  and  model  4  for  the  ac¬ 
celerating  flow.  Models  4  and  7  differ  in  the  same  way  as  models  2  and  6. 
For  the  afterbody  flow  an  incidental  comparison  is  made  in  Case  C  be¬ 
tween  models  4  and  3;  i.e.  ,  the  omission  of  the  +  O  chp.rge  transfer 
reaction  is  considered. 

Figure  10  shows  the  ionization  histories  along  the  stream  surface  passing 
through  the  stagnationpoint  for  Cases  A,  B,  and  C  relating  to  ballistic  re-entry 
conditions.  Case  A  exhibits  the  familiar  overshoot  of  electron  concentration 
causes  by  the  N  +  O  associative  ionization  reaction  seeking  a  local  equilibrium 
level  higher  than  the  final  level.  Characteristically,  there  is  little  change  from 
the  stagnation  region  to  the  sonic  point.  However,  the  dissociative  recombination 
is  rapid  enough  to  follow  the  rapidly  decreasing  local  equilibrium  in  the  supersonic 
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region  for  a  considerable  distance.  The  change  in  slope  in  the  shoulder  region 
is  caused  by  the  difference  between  the  decay  rates  of  the  Newtonian  and  blast 
wave  pressure  distributions  assumed.  It  was  brought  out  in  Ref.  19  that  the 
nonequilibrium  electron  concentrations  undershoot  the  infinite -rate  equilibrium 
curve. 

Case  B  exhibits  two  rather  distinct  regimes  in  the  shock  wave  relaxation 
zone.  The  early  portion,  at  low  ionization  levels,  is  controlled  by  impact  ioni¬ 
zation  and  the  final  rise,  by  associative  ionization.  It  should  be  observed  here 
that  the  rates  of  ionization  are  calculated  through  the  equilibrium  constant  by 
means  of  the  recombination  rate  constant  and  the  assumption  of  detailed  balancing. 

The  particular  reactions  governing  the  early  portion  of  the  rise  have  recombi- 

1 9 

nation  rate  constants  which  are  considerably  enlarged  by  inelasticity  factors. 

The  argument  for  using  such  factors  is  based  on  the  ability  of  molecular  colliders 
to  absorb  energy  in  internal  modes,  thereby  enhancing  greatly  the  probability 
of  recombination.  The  reverse  argument  for  the  ionization  rates  must,  there¬ 
fore,  be  based  on  the  possibility  of  a  vibrationally  excited  molecule  contributing 
its  internal  energy  to  help  ionize  a  collision  partner.  The  result  shown  in 
Fig.  10  is  based  on  the  assumption  of  vibrational  equilibration  throughout,  and 
must  be  regarded  as  an  upper  limit  for  the  effect  of  the  additional  reactions. 
Evidently,  this  is  not  as  bad  an  assumption  as  it  may  seem  because  fast  vibra¬ 
tional  excitation  paths  are  available  through  the  exchange  reactions  which  are 

37 

frequent  at  these  temperatures.  The  coupling  of  dissociation  with  vibration 

38 

as  discussed  by  Treanor  and  Marrone  is  not  significant  in  the  region  under  con¬ 
sideration  here. 

The  continuation  of  Case  B  through  the  overshoot  region  and  back  to  the 
sonic  line  does  not  show  significant  variance  with  the  qualitative  behavior  of 
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Case  A.  Levels  are  elevated  by  some  30%  because  of  the  inclusion  of  additional 
ions.  The  minimum  in  the  curve  for  Case  B  indicates  a  delay  in  the  ionization 
of  atomic  species  which  is  allowed  to  occur  only  by  impact  with  this  model. 
Beyond  the  sonic  line,  the  curve  for  Case  B  differs  markedly  from  that  for 
Case  A.  Since  atomic  ions  can  recombine  only  through  three-body  paths,  all  of 
the  atomic  ions  are  frozen  at  essentially  stagnation  levels.  This  is  typical  of 
ternary  reactions  for  the  flight  conditions  considered.  Thus,  when  all  of  the 
molecular  ions  are  neutralized,  the  remaining  electrons  freeze  at  a  concentra¬ 
tion  approximately  equal  to  the  sum  of  the  stagnation  levels  for  both  atomic  ions. 

Case  C  is  practically  identical  to  Case  D  throughout  the  decelerating  flow 
region.  The  charge  transfer  reactions  provide  a  more  efficient  means  of  pro¬ 
ducing  atomic  ions  downstream  from  the  overshoot  peak.  For  this  reason,  the 
dip  in  the  curve  preceding  the  stagnation  point  is  not  as  pronounced  as  it  is  with 
Case  B.  Downstream  from  the  sonic  line.  Case  C  more  nearly  resembles 
Case  A  than  it  does  Case  B.  The  atomic  ions  are  neutralized  by  receiving  elec¬ 
trons  from  neutral  molecules  through  the  charge  transfer  reactions.  The  re¬ 
action  chains  shown  in  Fig.  8  govern  the  system,  and  the  molecular  ions  re¬ 
combine  with  electrons  through  fast,  dissociative  reactions.  This  drives  the 
electron  concentration  below  that  in  Case  A  becan.se  the  local  equilibrium  ioni¬ 
zation  levels  are  actually  depressed  in  the  afterbody  region  by  the  inclusion  of 
the  additional  ions.  The  key  to  this  phenomena  lies  mainly  in  the  depression  of 
NO^  concentrations  at  stagnation  conditions.  Because  nitric  oxide  has  an  ioni¬ 
zation  potential  considerably  below  those  of  the  other  species  .and  because  there 
is  less  nitric  oxide  present,  the  local  equilibrium  ionization  curve  (which  is  es¬ 
sentially  determined  by  nonequiJibrium  behavior  of  the  neutral  species  which  are 
not  shown)  is  steeper  in  the  far  supersonic  portion  of  the  expanding  flow.  , 
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The  omission  of  reaction  33  which  cross-links  the  chains  shown  in  Fig.  8 

has  practically  no  effect  at  the  conditions  of  Fig.  10.  At  higher  altitudes  this 

+  -f 

omission  leads  to  sudden  freezing  of  N  because  the  N  +  N2  charge  transfer 
reaction  (28)  produces  at  a  greater  rate.  This  arises  from  higher  initial 
levels  of  N-atoms  and  earlier  freezing  of  N-atoms.  In  the  cases  studied,  initial 

N^- concentrations  are  higher  also.  The  sudden  freezing  then  is  caused  by  the 

« 

establishment  of  a  quasi- stationary  N^-condition  in  the  left-hand  branch  of  the 

charge  transfer  mechanism  shown  in  Fig.  8. 

It  has  been  demonstrated  by  these  comparisons  that  for  ballistic  re-entry 

conditions,  levels  and  gradients  of  electron  density  are  sensitive  to  the  choice 

of  reaction  mechanism  in  certain  regimes.  The  reflection  of  electromagnetic 
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waves  from  a  nonuniform  plasma  has  been  shown  ’  to  be  highly  dependent  upon 
both  these  features.  Therefore,  when  electromagnetic  propagation  problems 
are  considered,  it  would  appear  that  the  same  care  must  be  taken  in  selecting 
an  ionization  mechanism  as  has  been  taken  in  the  chemical  case.  The  impact 
ionization  reactions  may  play  an  important  role  in  the  early  part  of  the  shock 
wave  relaxation  zone,  and  the  charge  transfer  reactions  are  significant  in  the 
afterbody  nonequilibrium  flow.  Of  course,  the  associative  ionization  reactions 
and  their  reverse  reactions  are  basic  to  any  kinetic  scheme  for  ionized  species.. 
It  will  be  recalled  that  at  pressures  appropriate  to  shock  tunnel  operation,  the 
charge  transfer  reactions  eliminate  atomic  ions  early  in  the  expansion.  The 
remaining  portion  of  the  ionization  history  is,  therefore,  governed  by  the  simple 
model  1. 

3.  Lunar  Probe  Re-entry  Case 

The  ionization  history  for  a  lunar  probe  re-entry  flow  field  at  200  Kft  in 
altitude  is  dominated  by  atomic  ions.  Therefore,  the  complexities  of  model  7 
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must  be  included  until  some  rational  simplifications  are  discovered.  Figure  11 
shows  the  positive-ion  distributions  as  well  as  the  electron  distribution  for  such 
a  case.  The  molecular  ion  concentrations  display  large  overshoots  in  the  very 
early  part  of  the  shock  wave  relaxation  zone.  This  occurs  as  the  associative 
and  collisional  ionization  mechanisms  are  opposed  by  the  charge  transfer  re¬ 
actions  which  give  rise  to  the  atomic  ions.  Note  that  the  electron  concentration 
distribution  has  two  regimes;  one  corresponds  to  the  formation  of  molecular  ions 
and  the  other,  to  the  formation  of  atomic  ions.  This  explanation  differs  from 
that  given  for  the  ballistic  re-entry  case  which  is  dominated  throughout  by  mo¬ 
lecular  ion  formation  following  two  different  mechanisms  successively.  In  the 
afterbody  flow,  the  atomic  ions  freeze  leading  to  a  complete  freezing  of  ioniza¬ 
tion  at  a  relatively  higher  level  than  that  in  the  23  Kft/sep  case.  The  deioniza- 

« 

tion  chains  shown  in  Fig.  8  depend  on  the  presence  of  neutral  molecules  to  a 
great  extent.  If  concentrations  of  these  species  go  below  a  certain  level,  and 
molecular  ion  concentrations  are  low,  the  reaction  chains  do  not  function.  Note 
the  decrease  in  N^  concentration  and  the  increase  of  the  concentration  to 
the  final  levels.  This  occurs  because  of  the  different  behavior  of  the  local  equi¬ 
librium  trends  and  the  ultimate  freezing  of  the  charge  transfer  reactions. 

Other  information  relating  to  the  flow  field  is  given  in  Figs.  12,  13,  and 
14.  Neutral  particle  histories  shown  in  Fig.  12  qualitatively  resemble  those  ob¬ 
tained  for  lower  energy  re-entry  conditions.  Molecule  concentrations  are  gener¬ 
ally  much  lower,  except  nitric  oxide  overshoots  the  stagnation  equilibrium  level 
by  a  considerably  greater  amount.  Temperature  and  fluid  density  distributions 
are  shown  in  Figs,  13  and  14  respectively  since  these  are  quantities  of  interest 
in  plasma  and  radiation  calculations. 
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E.  Summary  of  Nonequilibriiim  Theoretical  Results 

The  establishment  of  further  generalization  in  ion  kinetic  models  has 
revealed  regimes  of  practical  interest  in  which  the  single  nitric  oxide  ionization 
model  may  not  be  sufficient  for  describing  the  electron  density  distributions  in 
nonequilibrium  airflows.  For  nozzle  flows  appropriate  to  high-enthalpy  aerody¬ 
namic  studies  at  reservoir  stagnation  temperatures  up  to  8000*K,  the  simple 
model  appears  to  be  adequate.  However,  in  the  early  part  of  the  relaxation  zone 
behind  the  strong  portions  of  a  detached  shock  wave  at  ballistic  re-entry 
conditions,  impact  ionization  of  molecules  may  be  significant.  In  the  afterbody 
and  wake  flows,  charge  transfer  processes  contribute  to  the  deionization 
mechanism.  For  superorbital  re-entry  flow,  a  complex  kinetic  model  may  be 
necessary  because  of  the  dominant  role  of  atomic  ions. 

The  present  work  must  he  rewarded  primarily  as  an  heuristic  enterprise 
at  the  present  time  because  of  important  unanswered  questions  in  high-temperature 
gas  kinetics.  Ionization  of  electronically  excited  species  and  excitation  kinetics 
have  not  been  included.  Vibrational- rotational  coupling  with  one  another  and 
with  dissociation  has  been  omitted.  Radiative  coupling  and  photo-processes 
have  not  been  considered  in  the  analysis.  Enhancement  of  three-body  processes 
due  to  electrons  as  catalytic  bodies  has  not  been  accounted  for.  Further 
experimental  research  guided  by  theory  must  be  undertaken  not  only  to  improve 
the  data  for  the  processes  included,  but  also  to  permit  the  rational,  considera¬ 
tion  of  the  additional  effects  enumerated  above. 
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IV.  IONIZATION  EXPERIMENTS  IN  THE  SHOCK  TUNNEL 


E 
¥ 


The  theoretical  results  given  in  the  previous  section  indicate  the  necessity 
of  rather  complex  ion-kinetic  models  for  certain  re-entry  body  flow  conditions. 
However,  they  further  show  that  model  I  (species  1-8,  reactions  1-10)  is  ade¬ 
quate  for  describing  nonequilibrium  ionization  in  shock  tunnel  flows  of  current 
interest.  In  effect,  the  electron  density  at  moderate-to-large  area  ratios  is 
determined  only  by  reaction  10  because  efficient  removal  processes  quickly  elim¬ 
inate  all  ions  except  NO^  and  e".  The  present  experimental  program  is  being 
pursued  in  an  effort  to  determine  the  rate  constant  for  reaction  10.  Previous 
experimental  studies  of  reaction  1  0  by  other  methods  have  been  discussed  in 
Section  III,  Part  B.  A  second  objective  of  the  experiment  was  the  study  of  the 
gasdynamic  processes  determining  the  duration  of  useful  test  flow  in  the  tunnel. 
Phase  shift  measurements  obtained  with  an  X-band  (9397  Me)  microwave  inter¬ 
ferometer  were  employed  to  obtain  electron  density  levels  at  the  2.  6-inch  en¬ 
trance  to  the  terminal  (conical)  stage  of  the  nozzle  system  (Fig.  2).  Pitot  pres¬ 
sure  measurements  were  obtained  simultaneously  in  the  six-foot  test  section  in 
an  effort  to  correlate  the  sequence  of  events  involved  in  starting  the  nozzle  with 
the  ionization  histories.  One  advantage  of  a  nozzle  kinetic  experiment  is  the 

removal  of  spatial  resolution  difficulties  which  may  be  encountered  in  shock 

,  20,21 

wave  relaxation  zone  experiments. 


A.  Microwave  Interferometer 

The  microwave  interferometer  was  used  to  measure  the  phase  shift  of 
1  electromagnetic  waves  due  to  propagation  through  the  plasma,  which  acts  as  a 

r  dielectric  medium.  This  technique  has  been  widely  applied  in  gaseous  discharge 

;  research  and  has  been  described  fully  in  the  literature  relating  to  that  field 

r 
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(e.g.  Refs.  39  and  40).  Microwave  diagnostic  methods  are  analogous  to  optical 
probes  in  that  the  propagating  signal  is  so  weak  that  it  leaves  the  test  gas  es¬ 
sentially  undisturbed. 

Figure  15  shows  a  schematic  diagram  of  the  bridge  circuit  of  which  the 
interferometer  is  comprised.  Power  is  supplied  by  a  klystron  oscillator  through 
an  isolator  in  order  to  prevent  reflected  power  from  entering  the  oscillator,  A 
directional  coupler  is  employed  to  divide  the  po'^er  between  the  test  arm  and  the 
reference  arm.  Cross-guide  couplers  with  crystal  detectors  are  mounted  on 
either  side  of  the  test  region  in  order  to  measure  reflected  and  transmitted 
power.  Horns  are  used  to  obtain  transmission  through  the  plasma.  The  refer¬ 
ence  arm  is  equipped  with  a  phase-shifter  and  an  attenuator  which  are  used  both 
to  initially  balance  the  bridge  and  to  quantitatively  measure  the  unbalance  pro¬ 
duced  by  the  plasma.  A  detector  is  used  to  monitor  power  level  in  the  reference 
arm.  A  magic-tee  is  used  to  mix  the  power  from  the  two  .-jrms.  The  mixed 

signal  is  fed  to  the  final  detector  through  an  isolator.  Thus,  the  detected  mixed 

* 

signal  is  affected  by  attenuation  and  phase- shift  in  the  plasma. 

Figure  16  shows  the  physical  arrangement  and  location  of  the  microwave 
circuitry.  The  transmitted  signal  is  directed  normal  to  flow  direction  through 
windows  made  of  bonded  glass  fiber.  The  thickness  of  the  windows  was  selected 
for  optimal  transmission  of  X-band  (9397  Me)  microwave  radiation.  The  test 
signal  passes  through  the  plasma  in  a  region  ynmediately  downstream  from  an 
expansion  fan  which  turns  the  flow  through  10'  and  upstream  from  the  entrance 
lips  to  the  conical  nozzle.  Turning  the  flow  centrifuges  out  any  solid  material. 
Therefore,  the  flow  is  gasdynamically  uniform  with  the  exception  of  the  boundary 
layers  on  the  windows.  Boundary -layer  growth  in  the  primary  nozzle  was  com¬ 
puted  (for  the  range  of  conditions  used)  by  means  of  an  empirical  relationship 
given  in  Ref.  41,  The  displacement  thickness  of  the  boundary  layer  on  the 
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Fiberglas  windows  was  approximately  0.4  of  a  free-space  wavelength  at  75  atm 
reservoir  pressure  and  decreased  to  0.15  at  680  atm  while  the  total  path  length 
through  the  plasma  is  about  four  free-space  wavelengths  at  X-band. 

The  first  step  in  calibrating  the  interferometer  is  tuning  out  the  mismatch 
between  the  horns  with  no  plasma  present.  This  is  necessary  because  multiple 
reflections  between  the  horns  give  rise  to  a  standing  wave.  It  is  difficult  to 
maintain  near-perfect  tuning  when  a  weak  plasma  is  introduced  because  small 
alterations  in  the  propagation  path  change  the  strength  of  the  horns  standing 
wave,  and  this  produces  either  a  loss  or  even  a  gain  in  the  transmitted  signal. 
Before  a  test,  the  phase -amplitude  bridge  is  nulled  by  adjusting  the  reference 
phase  shift  and  attenuation  until  no  signal  is  detected  by  the  mixed  signal  de¬ 
tecting  crystal.  In  order  to  test  the  tuning  and  the  interaction  between  the  two 
arms,  carefully  fitted  polymer  foam  blocks  having  known  dielectric  and  loss 
properties  were  placed  between  the  windows  in  the  nozzle.  Measured  phase 
shifts  agreed  with  predicted  values  to  within  10%.  The  shock  tube  supplying  the 
tunnel  was  fitted  with  a  closure  plate  which  blanked  off  the  nozzle  entrance  and 
was  operated  at  conditions  typical  of  those  used  for  nozzle  runs.  Interfero¬ 
metric  data  were  then  obtained  to  determine  the  disturbance  levels  of  recoil 
motion  and  of  elastic  wave  propagation.  It  was  found  tha  t  these  mechanical 
effects  could  be  safely  neglected  in  the  interpretation  of  the  data  , 

B .  Nonequilibrium  Ionization  Measurements 

The  shock  tube  supplying  the  tunnel  was  operated  with  hydrogen  driver  gas 
at  pressures  from  500  psi  to  12,  000  psi  and  over  an  incident  shock  Mach  number 
range  of  6,  7  to  10.  7.  These  conditions  gave  a  range  of  reservoir  states  from 
T  =  3500°K,  p  =  25  atm  to  T  =  7100'’K,  p  =  650  atm.  The  low  pressures 

O  O  0  ^0 

and  Mach  numbers  were  selected  to  minimize  the  air  ionization  at  the  test 
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station.  Therefore  the  electron  density  levels  which  were  measured  (Fig.  17) 

serve  as  a  measure  of  the  maximum  impurity  ionization  level.  Ionization  curves 

1 3 

for  sodium  (I  =  5. 12  ev)  show  that  fo-^r  10  atoms/cc  70%  of  the  atoms  are  io¬ 
nized  at  3500*K  and  that  there  is  complete  ionization  at  4000‘‘K .  Conseq\iently, 
at  all  higher  temperatures  there  should  be  no  increase  in  the  number  of  im¬ 
purity  electrons. 

* 

Microwave  records  were  obtained  for  reference  arm  power,  plasma  re¬ 
flected  power,  plasma  transmitted  power,  and  mixed  signal  power.  The  signal 
level  through  the  reference  arm  remained  constant  during  all  runs  indicating  no 
unwanted  interaction  of  this  signal  with  the  plasma  in  the  test  arm.  The  power 
reflection  coefficient  of  the  plasma  was  less  than  0. 1%  for  all  experiments.  The 
data  were  reduced  to  phase- shift  readings  by  adjusting  the  phase- shifter  and  at¬ 
tenuator  in  the  reference  arm  to  duplicate  the  combination  of  transmitted  power 
level  and  mixed  signal  level  recorded  at  each  time  during  the  experiment.  At¬ 
tenuation  due  to  the  plasma  could  be  obtained  because  it  was  larger  than  the  horn 
detuning  effect  described  above.  Some  transmission  loss  may  be  accountable 
on  the  basis  of  near-field  radiation  pattern  distortion  due  to  slight  transverse 
nonuniformity  of  the  plasma.  The  magnitude  of  this  gradient,  which  is  perpen¬ 
dicular  both  to  the  flow  and  the  microwave  beam,  has  been  calculated  to  give  no 
more  than  4%  variation  about  the  mean  electron  density.  In  the  following  sub¬ 
section  C  on  theoretical  interpretation,  both  the  origin  and  the  method  of  pre¬ 
diction  of  this  nonuniformity  are  discussed  in  greater  detail.  Measurement  of 
attenuation  is  more  sensitive  to  this  type  of  distortion  than  measurement  of 
phase  shift  because  of  the  directionality  of  the  horns.  However,  auxiliary  de¬ 
tectors  situated  at  many  locations  failed  to  detect  any  refracted  signal  during 
the  runs. 
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Each  phase- shift  record  obtained  from  the  reduction  procedure  described 
above  first  exhibited  a  sharp  rise  signifying  the  arrival  of  the  leading  compres¬ 
sion  wave  preceding  the  gas  originating  in  the  reservoir.  The  rise  was  followed 
by  an  irregular  but  high  level  of  signal  which  dropped  to  a  plateau  having  a  com¬ 
paratively  low  noise  level  when  the  (u-a)  wave  arrived  signaling  the  end  of  the 
tunnel  starting.  This  plateau  signified  the  useful,  quasi-steady  test  flow  of  air. 
After  the  first  plateau  the  ionization  dropped  to  a  second  plateau  at  a  level  about 
one-half  that  of  the  first.  Finally,  the  level  dropped  back  to  zero.  The  end  of 
the  useful  test  flow  is  marked  by  the  first  measurable  effect  of  driver  gas  di¬ 
lution  of  the  air  terminating  the  first  plateau.  It  seems  plausible  that  the 
second  plateau  arises  from  the  nonequilibrium  ionization  frozen  in  the  quenched 
combustion  zone  originating  from  an  air -hydrogen  diffusion  flame  sustained  by 
interfacial  mixing, 

C.  Interpretation  of  the  Results 

1 .  Establishment  of  Useful  Testing  Times 

The  explanation  of  events  in  the  ionization  histories  was  qualitatively  sup¬ 
ported  by  theoretical  estimates  of  the  timing  of  the  starting  flow.  The  starting 
process  in  a  shock  tunnel  operating  behind  the  reflected  shock  wave  is  discussed 
in  Ref.  42.  That  report  shows  that  the  tunnel  will  start  rapidly  when  the  initial 
pressure  in  the  nozzle  is  very  low.  The  ideal  minimum  time  for  steady  flow  to 
be  established  is  then  the  time  required  for  a  nonsteady  expansion  fan  which  is 
initiated  at  the  throat  to  sweep  pass  the  testing  station.  Combining  this  compu¬ 
ted  time  with  the  time  required  for  the  incident  shock  wave  to  travel  from  the 
triggering  station  to  the  end  wall  consistently  predicted  a  starting  time  for 
steady  flow  which  was  approximately  100  fisec  less  than  the  measured  value. 

This  discrepancy  is  not  surprising  because  the  actual  starting  process  is  two- 
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dimensional  with  waves  propagating  from  the  walls  through  the  boundary  layer 
to  the  center  of  the  flow  as  well  as  propagating  axially  along  the  nozzle.  This 
difference  in  starting  time  became  very  noticeable  when  a  correlation  with  heat 
transfer  and  stagnation  pressure  in  the  six  foot  test  section  was  made.  At  this 
station  the  time  difference  was  the  order  of  several  hundred  microseconds. 

Uncertainties  due  to  interface  mixing  prevent  a  prediction  of  the  end  of 
42 

useful  test  time.  However,  the  ionization  records  showed  rather  sudden 
changes  to  negative  slopes  at  the  end  of  the  first  plateau.  The  interface  travels 
down  the  nozzle  on  a  u-trajectory  proceeded  by  a  (u  +  a)  wave;  however,  because 
of  the  low  temperatures,  both  the  wave  and  interface  arrive  in  the  test  section 
at  nearly  the  same  time.  Extrapolation  of  the  end  of  steady  flow  as  measured 
by  the  microwave  interferometer  to  the  six  foot  test  section  again  indicated  a 
discrepancy  with  the  heat  transfer  and  stagnation  pressure  of  several  hundred 
microseconds. 

While  neither  the  correlation  of  the  beginning  or  end  of  the  useful  test 
time  in  the  test  section  can  be  made  on  the  basis  of  a  simple  one-dimensional 
wave  analysis,  the  agreement  of  the  total  time  difference  from  beginning  to  end 
was  in  general  fairly  good.  The  actual  testing  times  were  much  less  than  that 
predicted  on  the  basis  of  an  ideal  interface.  However,  useful  test  flow  was 
achieved  over  a  wide  range  of  conditions.  The  trend  of  the  data  indicates  that 
testing  time  increases  somewhat  as  the  reservoir  pressure  level  increases  and 
that  it  decreases  as  the  Mach  number  increases  for  Mach  numbers  between  6.7 
and  10.7.  Extrapolation  of  test  section  testing  time  versus  incident  shock 
strength  to  zero  time  gives  an  intercept  at  a  Mach  number  in  the  vicinity  of  11. 
These  results  indicate  that  the  interface  which  terminates  the  steady-flow 
region  arrives  far  earlier  than  ideal  shock-tube  theory  predicts  (based  on 
planar  interface). 
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2.  Interface  Mixing  Phenomena 

As  pointed  out  in  the  preceding  section  the  availability  of  the  microwave 
interferometer  technique  has  proved  a  reliable  method  of  determining  the 
arrival  of  the  interface  in  the  six-foot  shock  tunnel.  Previously,  reliance  was 
placed  on  a  combination  of  pitot  static  and  heat  transfer  instrumentation  to 
establish  interface  position.  With  high-enthalpy  drivers  this  has  not  always 
proved  reliable.  With  microwave  instrumentation  the  observation  of  the  depar¬ 
ture  of  available  testing  time  in  the  six-foot  shock  tunnel  proved  to  be  at  variance 
with  previous  experience  obtained  in  shock  tunnels  at  Cornell  Aeronautical 
Laboratory.  Previously,  utilizing  cold  helium  as  the  driver  in  the  11"  by  15" 
shock  tunnel,  better  than  60%  of  the  available  test  time  could  be  achieved,  and 
even  with  cold  hydrogen  at  least  50%  of  the  test  time  would  be  achieved.  ^  In 
general,  the  pressures  within  the  shock  tube  were  very  high  during  shock  tunnel 

operation  and  this  reduction  in  testing  time  was  attributed  to  mixing  initiated 

43 

by  rupturing  of  the  diaphragm  rather  than  to  viscous  phenomena.  The  reduc¬ 
tion  of  testing  time  in  the  present  shock  tunnel  is  a  result  of  complex  gasdynamic 
phenomena  and  may  be  directly  related  to  the  specific  design  of  the  six-foot 
shock  tunnel.  In  the  succeeding  paragraphs,  an  explanation  is  offered:  however, 
until  definitive  experiments  are  carried  out  this  should  be  regarded  as  tentative. 

As  indicated  in  the  description  of  the  shock  tunnel  facility,  heated  hydro¬ 
gen  is  utilized  to  obtain  tailoring  for  high-enthalpy  operation.  The  performance 
indicated  that  ideal  shock-tube  theory  must  be  modified  when  utilizing  heated 
drivers.  In  this  shock  tube,  the  driver  is  heated  externally  prior  to  a  run 
in  order  to  avoid  the  development  of  internal  thermal  stresses.  In  order  to 
prevent  damage  to  the  O- rings  sealing  the  diaphragm,  cooling  coils  are 
placed  adjacent  to  the  diaphragm  region.  This  leads  to  a  temperature 
gradient  within  the  driver  so  that  prior  to  each  run,  cold  hydrogen 
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is  collected  in  the  vicinity  of  the  diaphragm.  Immediately  after  bursting  of  the 
diaphragm,  only  cold  hydrogen  is  available  as  the  driver  gas.  Subsequently, 
as  the  cold  hydrogen  is  used  up,  it  must  in  turn  be  accelerated  by  the  heated 
hydrogen.  This  difficulty  is  further  enhanced  by  the  large  degree  of  chambrage 
{area  contraction)  incorporated  into  the  shock  tube  design  with  the  hope  of  in¬ 
creasing  performance.  The  area  contraction  lowers  the  velocity  within  the 
driver,  and  a  longer  time  is  required  before  the  cold  hydrogen  is  swept  out  of 
the  driver.  Thus,  the  driven  gas  is  effectively  accelerated  by  a  driver  gas  with 
an  axial  temperature  gradient. 

This  type  d  r  i  ve  r  introduces  new  and  unanticipated  problems  in 
relationship  to  tailoring.  Tailoring  requires  that  the  driver-driven  interface 
be  so  matched  that  the  reflection  of  the  shock  wave  from  the  nozzle  end  does 
not  create  further  reflections.  In  the  case  of  a  driver  of  uniform  temperature, 
the  tailoring  conditions  can  be  directly  calculated.  However,  with  the  tempera¬ 
ture  gradient  in  the  driver,  no  formal  tailoring  condition  can  be  found.  The 

interface  between  the  air  and  the  hydrogen  generally  involves  a  density  gradient 

44 

which  is  expected  to  be  unfavorable  from  the  viewpoint  of  Taylor  stability, 

and  this  condition  is  further  complicated  by  the  observation  of  combustion  at 

the  interface.  Tailoring,  utilizing  the  heated  hydrogen  driver,  is  not  normally 

calculated  but  is  based  on  observation  of  the  pressure  history  in  the  shock  tube 

and  the  examination  of  the  available  testing  time.  Also,  equilibrium  interface 
7  45 

operation  ’  is  often  used.  This  technique  employs  secondary  reflections 
from  the  interface  to  further  increase  the  enthalpy,  and  this  operating  condition 
tends  to  produce  regions  which  promote  Taylor  instability.  Careful  monitoring 
of  the  flow  then  becomes  mandatory. 

In  view  of  the  complexity  of  the  phenomena  assi'ciaterl  with  the  reduction 
of  testing  time  during  current  operations  of  this  shock  tunnel,  it  would  be 
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unwise  at  the  present  time  to  generalize  about  the  factors  involved  in  this  testing 
time  loss.  Currently  a  research  program  is  underway  to  study  the  aspects  of 
tailored  instability  as  they  affect  interface  mixing,  as  well  as  a  program  of 
modification  to  the  shock  tunnel  in  order  to  extend  its  range  of  operation. 


3.  Deionization  Rate  Constant  Experiment 

19  33 

As  demonstrated  previously,  ’  the  electron  density  at  a  given  area 
ratio  in  nonequilibrium  expanding  airflows  is  relatively  insensitive  to  reservoir 
pressure,  but  is  almost  inversely  proportional  to  the  rate  constant.  Estimates 
of  the  electron  density  were  made  with  the  aid  of  numerical  solutions  of  non¬ 
equilibrium  nozzle  flows  from  reservoir  temperatures  from  5000‘’K  to  8250®K 
at  a  few  typical  reservoir  pressures.  The  calculations  were  performed  using 
the  program  described  in  Ref.  32  and  the  neutral  particle  kinetic  data  shown 
in  Table  1.  Reaction  10  (forward  and  reverse)  shown  in  Table  2  was  the  only 
ionization  reaction  considered  for  the  reasons  given  at  the  beginning  of  this 
section.  A  polynomial  curve  fit  of  the  geometrical  area  function  of  the  primary 
nozzle  contour  was  employed.  Boundary-layer  corrections  based  on  the  methods 
of  Ref.  41  were  applied  to  get  the  gasdynamic  conditions  at  the  exit  of  the 
nozzle;  however,  the  electron  mass  concentration  at  the  geometrical  exit  area 
was  used.  This  is  a  better  approximation  than  the  effective  area  ratio  because 
the  slope  of  the  area  function  is  small  near  the  exit.  Hence,  the  integrated 
residence  time  would  be  more  nearly  described  by  choosing  the  nonequilibrium 
ionization  level  based  on  actual  axial  distance  rather  than  effective  area  ratio. 

Since  the  chemistry  of  the  neutral  particles  was  essentially  frozen  at  the 
conditions  of  the  experiment,  ideal  gas  calcvilations  were  employed  for  ob¬ 
taining  conditions  after  the  expansion  fan  which  turns  the  flow  through  10°. 
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Frozen  vibration  was  assumed  to  obtain  the  effective  specific  heat  ratio.  Ioni¬ 
zation  relaxation  was  calculated  for  the  inner  and  outer  streamlines  entering 
the  microwave  beam.  Line  segment  approximations  with  average  densities 
were  used  for  the  outer  streamline.  Computation  of  ionization  relaxation  for 
the  inner  streamline  was  straightforward  since  the  residence  time  in  the  expan¬ 
sion  fan  is  approaching  zero.  Over-all  reductions  in  electron  mass  concentra¬ 
tions  between  primary  nozzle  exit  and  the  microwave  beam  were  in  the  near 
vicinity  of  30%  for  all  cases.  Greater  recombination  along  the  outer  stream¬ 
line  than  the  inner  streamline  resulted  in  a  plasma  gradient  which  was  mutually 
perpendicular  to  the  flow  velocity  and  the  microwave  propagation  direction.  As 
mentioned  previously,  this  variation  never  resulted  in  nonuniformity  exceeding 
4%  about  the  mean  value.  The  extent  of  recombination  in  the  expansion  fan  was 
rather  insensitive  to  the  computational  model,  the refore,  no  higher  approxima¬ 
tions  for  the  outer  streamline  were  sought. 

Figures  17  and  18  show  the  preliminary  data  obtained  from  llie  experi¬ 
ments.  The  error  bars  are  based  only  on  the  uncertainty  in  reading  the  oscil¬ 
loscope  traces.  They  do  not  include  effects  of  detuning,  error,s  in  the  calibra¬ 
tion  procedure,  assumption  of  a  one  dimen-sional  microwave  In'ini  for  data  re¬ 
duction,  or  neglect  of  electron  density  gradients  in  tlie  boundary  layer.  The 
electron  density  was  found  to  be  independent  of  the  pressure  level.  This  in¬ 
sensitivity  was  also  exhibited  by  the  results  from  the  complete  lumequilibr turn 
flow  solutions. 

The  atom  and  molecule  cross-sections  used  for  predicting  the  electron 
collision  frequencies  in  Fig,  18  were  taken  from  F'ig  7  in  Ref.  46.  The  mea- 
svired  electron  collision  frequencies  tend  to  be  higher  than  those  predicted. 

This  difference  was  most  pronounced  for  the  runs  where  the  pressure  varied 
between  75  and  100  atmospheres.  Assuming  that  the  aver.igi'  temper, iture  of 
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the  electrons  corresponds  to  the  local  gas  temperature  the  collision  frequency 

g 

•would  be  4  X  10  cps,  however,  if  the  electron  temperature  was  at  the  temper¬ 
ature  the  gas  had  at  the  nozzle  throat  (an  approximation  to  the  gas  vibrational 

temperature  at  the  testing  station)  the  collision  frequency  would  have  been 
9  9 

1 . 3  X  10  cps.  The  measured  value  of  3  x  10  cps  fell  outside  of  this  range. 

Corrections  to  the  collision  frequency  which  account  for  electron-ion  collisions, 

7  4  13 

2x10  cps,  electron- sodium  atom  collisions,  10  cps  (assviming  10  atoms/ 

7 

cc  in  the  reservoir),  and  electron- copper  atom  collisions,  10  cps  (assuming 

0.1  cc  of  nozzle  throat  is  vaporized,  an  extreme  uppei;  limit)  are  incapable  of 

accounting  for  the  differences.  The  effect  of  an  electron  density  gradients  in 

the  boundary  layer  would  be  to  decrease  both  the  phase  change  and  the  attenu- 

47 

ation  of  the  transmitted  signal.  Therefore  the  measured  electron  collision 
frequency  in  Fig.  18  would  be  decreased  and  the  electron  densities  in  Fig.  17 
would,  in  general,  be  increased  if  this  effect  were  taken  into  piccount.  The 
quantitative  effect  of  this  gradient  on  the  reduced  data  is  at  present  not  known. 
The  main  difficulty  in  correcting  the  data  is  in  determining  the  distribution  of 
electrons  in  the  boundary  layer.  The  more  direct  method  would  be  to  employ 
boundary  layer  bleeds  just  upstream  of  the  measuring  station  and  thus  to  avoid 
making  this  an  important  correction  to  the  data. 

The  value  of  the  dissociative-recombination  rate  constant  which  gives 
good  agreement  with  the  measured  electron  densities  is  of  the  same  order  of 
magnitude  as  that  in  Ref.  18.  A  more  precise  determination  of  this  rate  con¬ 
stant  will  be  made  following  additional  experiments  which  are  now  underway. 
Specifically,  runs  with  nitrogen  instead  of  air  as  the  test  gas  are  being  made  to 
verify  that  the  measured  phase  shift  and  attenuations  are  not  due  to  impurities. 
Furthermore  a  simultaneous  reading  S-band  interferometer  placed  near  the 
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entrance  of  the  Fiberglas  portion  of  the  conical  nozzle  should  yield  further  in¬ 
formation  on  the  temperature  dependence  of  the  rate  constant.  Extension  of  the 
theoretical  studies  may  shed  additional  light  on  the  questions  of  unimolecular 
processes  posed  in  Ref.  21  because  of  the  higher  density  levels  in  the  present 
experiment.  Values  of  the  specific  rate  constant  will  be  reported  after  further 
studies  are  made. 

In  summary,  the  microwave  diagnostic  methods  commonly  used  in  gaseous 
discharge  physics  have  proved  to  be  useful  tools  in  studying  not  only  nonequili¬ 
brium  ionization  in  shock  tunnel  flows,  but  also  the  basic  gasdynamic  processes, 
occuring  in  such  flows.  Considering  the  state  of  knowledge  in  high-temperature 

chemical  kinetics,  the  rate  information  indicated  by  the  preliminary  data  does 

20  21  48 

not  disagree  seriously  with  that  reported  elsewhere.  ’  ’ 
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V.  SUMMARIZING  REMARKS 


Shock  tunnel  methods  have  been  discussed  in  the  context  of  ionized  flow 
studies.  Certain  features  of  the  six-foot  facility  at  the  Cornell  Aeronautical 
Laboratory  have  been  described.  The  large  test  section  and  the  Fiberglas 
nozzle  render  it  useful  in  the  investigation  of  plasma  flows  appropriate  to  hyper¬ 
sonic  flight. 

One  of  the  most  significant  considerations  in  such  shock  tunnel  research 
is  that  of  nonequilibrium  phenomena  both  in  the  tunnel  flow  and  in  the  flow  near 
a  body.  Emphasis  in  the  present  work  has  been  on  ionization  mechanisms  in 
these  nonequilibrium  flows.  For  certain  purposes  of  electromagnetic  studies, 
rather  complex  mechanisms  may  be  necessary  for  blunt-body  flows  of  present 
day  engineering  interest  occurring  under  ballistic  re-entry  conditions.  Super- 
orbital  flight  poses  re-entry  problems  involving  large  amounts  of  energy  in¬ 
vested  in  ionization  and  requiring  the  use  of  generalized  ion-kinetics  for  real¬ 
istic  determination  of  flow-field  gasdynamics.  At  pressures  and  temperatures 
encountered  in  nozzle  flows  in  existing  shock  tunnel  facilities,  the  compara¬ 
tively  simple  ionization  models  presently  in  use  seem  to  suffice  for  large 
expansions. 

Microwave  interferometric  measurements  of  electron  density  have  been 
described.  The  records  indicate  that  the  effects  of  the  driver  gas  terminate 
the  useful  test  flow  earlier  than  would  be  predicted  on  the  basis  of  simple  shock 
tube  theory.  The  simplicity  of  ion  kinetics  in  the  shock  tunnel  flow  permits  the 
study  of  a  single  deionization  process.  Preliminary  data  agree  reasonably  well 
with  previously  reported  values  of  the  dissociative  recombination  rate  constant 
for  nitric  oxide  ions  and  electrons,  however  the  quantitative  effects  of  impuri¬ 
ties  in  this  facility  have  not  yet  been  fully  determined. 
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SPECIES  NUMBER  KEY  *N  s  ANY  COLLISION  PARTNER 


TABLE  HI 

CHARGE  TRANSFER  RATE  CONSTANTS 
COMPUTED  FROM  THEORY  OF  REFERENCE  23 


REACTION 

i 

REDUCED  HASS 

jU 

amu 

POLARIZABILITY 
cc  X  10*® 

cm® 

RATE  CONSTANT 

l<F|Xlo'° 

cm®  PARTICLE*'  sec*' 

23 

■raEH 

80 

6.68 

24 

176 

9.75 

2S 

160 

9.10 

28 

128 

8.66 

29 

160 

9.51 

30 

9.74 

160 

9.51 

33 

7.47 

80 

7.64 

TABLE  H 

RATIOS  OF  ACTUAL  TO  EQUILIBRIUM  AIRFLOW  VARIABLES  FOR 
COUPLED-REACTION  EXPANSIONS  GIVING  EQUILIBRIUM  MACH  NUMBER  20 

REF.  32 

^  =  I  cm 


To 

6000'’K 

«000''K 

"o 

100  atm 

1000  atm 

100  atm 

1000  etm 

0.425 

0.875 

0.159 

0.544 

T/Teq 

0.359 

0.855 

0.121 

0.497 

a^eq 

ITPM 

1.014 

I.IIO 

1.045 

"A  EQ 

0.913 

0.987 

0.901 

mm 

M/Meq 

1.455 

1.06 

2.30 

1.305 

Re, /Re, 

X'  Xgp 

2.79 

1.17 

8.25 

2.01 
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Pigurt  5  TEST  SEaiON  PITOT  PRESSURE  DISTRIBUTION 
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Figure  8  DEIONIZaTiON  REACTION  CHAINS  FOR  ELECTRONS  AND  SINGLY 
CHARGED  rOSrriVE-IONS  IN  AN  AIR  PLASMA 


(innixm  JO  II4VXO  mu  IllOV^t)  NOliVViNlONOO  NOMiOaii 
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Figure  9  ELECTRON  CONCENTRATION  DISTRIBUTIONS  FOR  AIRFLOW  IN  A  HYPERSONIC 

T_  =  8000”K;  ^  100  ofm;  -£  =  1  cm 
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Figun  10  lONIZATSON  PROFILES  ALONG  A  STREAM  SURFACE  ABOUT  A  HEMISPHERE-CYLINDER 

VEL.  =  23  kft/sec;  ALT.  =  150  kft;  IL  =  lit  (A,B  &  C  REFER  TO  CASES  DEFINED  IN  THE  TEXT) 


STAGMATtON!  SONIC 
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0  .01  .02  .03  .04  .05 

Figure  13a  TEMPERATURE  DISTRIBUTION  ALONG  STAGNATION 
STREAMLINE  IN  THE  SHOCK  WAVE  DETACHMENT 
ZONE  OF  A  HEMISPHERE-CYLINDER  BODY  (VEL  »  36 
Kfv^te;  ALT  200  Kit;  Pg  1  foot,  ZL  »  249»K) 
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Pigurt  14a  DENSITY  DISTRIBUTION  ALONG  STAGNATION  STREAMLINE  IN 
THE  SHOCK  WAVE  DETACHMENT  ZONE  OF  A  HEMISPHERE 
CYLINDER  BODY  (VEL  «  36  Kft/see;  ALT  =  200  Kft; 

1  foot),  =  3.18 X  lO-^gm/em^) 
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Figure  14b  DENSITY  DiSTRiBUTICN  ABOUT  NOSE  AND  AFTERBODY  OF  A 
HEMISPHERE.CYLINDER  (VEL  =  36  Kff/sec;  ALT  »  20C  Kft; 
^*0  *  1  foot;  ■  3.18  X  10*^  gm/cm^) 


Figure  15a  MICROWAVE  INTERFEROMETER  SCHEMATIC 

!| 
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Figure  15b  MICROWAVE  INTERFEROMETER  LOCATION 
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RESERVOIR  TEMPERRTURE.  ‘‘K 

Figure  18  ELECTRON  COLLISION  FREQUENCY  AT  MICROWAVE  STATION 
IN  SHOCK  TUNNEL  NOZZLE 


